Interplay between static and dynamic polar correlations in relaxor Pb(Mgi/ 3 Nb 2 /3)03 
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We have characterized the dynamics of the polar nanoregions in Pb(Mg 1/ / 3 Nb 2 /3)03 (PMN) 
through high-resolution neutron backscattering and spin-echo measurements of the diffuse scat- 
tering cross section. We find that the diffuse scattering intensity consists of both static and dynamic 
components. The static component first appears at the Curie temperature O ~ 400 K, while the 
dynamic component freezes completely at the temperature T/ ~ 200 K; together, these compo- 
nents account for all of the observed spectral weight contributing to the diffuse scattering cross 
section. The integrated intensity of the dynamic component peaks near the temperature at which 
the frequency-dependent dielectric constant reaches a maximum (T m ax) when measured at 1 GHz, 
i. e. on a timescale of ~ 1 ns. Our neutron scattering results can thus be directly related to dielectric 
and infra-red measurements of the polar nanoregions. Finally, the global temperature dependence of 
the diffuse scattering can be understood in terms of just two temperature scales, which is consistent 
with random field models. 

PACS numbers: 74.72.-h, 75.25,+z, 75.40.Gb 



I. INTRODUCTION 



Lead-based relaxor ferroelectrics are compounds with 
the general formula PbB0 3 in which the B-site is dis- 
ordered. They exhibit exceptionally large piezoelectric 
coefficients (d 3 3 ~ 2,500pC/N) and dielectric constants 
(e ~ 25, 000) making them attractive for device applica- 
tions^-. 3 . PbMg 1/3 Nb 2/3 3 (PMN) and PbZn 1/3 Nb 2/3 03 
(PZN) are prototypical relaxor ferroelectrics and the 
most studied; both display a broad and unusually 
frequency-dependent zero-field dielectric response (see 
Fig. [TJ>) that contrasts with the sharp and (compara- 
tively) frequency-independent peaks observed in conven- 
tional ferroelectrics such as PbTi0 3 (PT) and BaTi0 3 . 

This anomalous dielectric behavior is matched by the 
odd structural properties of the lead-based relaxors^ In- 
stead of a well-defined structural transition to a long- 
range ordered ferroelectric ground state, which normally 
characterizes a typical ferroelectric, lead-based relaxors 
develop short-range ferroelectric correlations on cooling 
that are consistent with tiny domains of ferroelectric or- 
der embedded within a paraelectric matrix, while the av- 
erage structural unit cell remains cubic. These local re- 
gions of ferroelectric order, now known as polar nanore- 
gions or PNR, were first postulated to explain the tem- 
perature dependence of the optical index of refraction 
of a variety of disordered ferroelectric materials j 5 - The 
existence of PNR has since been confirmed by numer- 
ous x-ray and neutron scattering studies, which report 
the presence of strong diffuse scattering at low temper- 
atures*^^ The diffuse scattering in PMN, for exam- 



ple, was recently investigated in great detail and shown 
to vanish above ~ 420 KJZ Typical diffuse scattering 
intensity contours measured at 200 K are displayed in 
Fig. [2Jz)i2 and show that the neutron diffuse scattering 
cross section is very broad in momentum, which implies 
the existence of short-range structural correlations. By 
comparison, the sharp, resolution-limited Bragg peaks 
that accompany a transition to a normal ferroelectric 
ground state such as in PbTi0 3 are indicative of long- 
range structural correlations. While the lead-based re- 
laxor diffuse scattering cross section has been modeled 
extensively, considerable debate persists over the under- 
lying physical origin of the intriguing butterfly-shaped 
contours illustrated in Fig.r2Jz)^^— 

Despite uncertainties in the origin of the diffuse scat- 
tering cross section, the polar nature of the diffuse scat- 
tering in the lead-based relaxors has been well estab- 
lished through three different means. First, the onset 
of the diffuse scattering in PMN was shown to coincide 
with the temperature at which a ferroelectric-active, soft, 
transverse optic mode reaches a minimum frequency (see 
Fig-IDO > an d which also coincides with the Curie- Weiss 
temperature derived from high-temperature susceptibil- 
ity measurements (Fig.[T]c)ii£~— Second, the diffuse scat- 
tering was shown to respond strongly to an electric field, 
which in PMN suppresses the diffuse scattering while si- 
multaneously enhancing the Bragg peaks at low temper- 
aturesi 2 ^ 3 - Third, on doping PMN with PT both the 
piezoelectric properties and the total diffuse scattering 
increase but then drop sharply above PT concentrations 
near 32%, at which point the diffuse scattering is replaced 
by a well-defined structural transition. 24 ' 25 While there is 
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a clear connection between the anomalous dielectric prop- 
erties and the diffuse scattering in lead-based relaxors, 
the lattice dynamics remain poorly understood, particu- 
larly at long timescales. Specifically, broadband infrared 
and frequency-dependent dielectric measurements, and 
the diffuse scattering cross section measured with neu- 
trons and x-rays are difficult to reconcile because the for- 
mer two techniques generally probe only the momentum 
response near Q = 0. This problem is further compli- 
cated by the fact that neutron scattering measurements 
are typically limited to energy resolutions SE ~ ITHz, 
and therefore direct comparisons with low-frequency di- 
electric data have not been possible. 

While comparisons of dielectric and infra-red measure- 
ments to neutron and x-ray scattering data are difficult 
to make, data at the high frequency limit of broadband 
measurements (~ 1 THz) on thin films have proven to 
be in excellent agreement with neutron inelastic scatter- 
ing measurements of the soft transverse optic mode. 26 
Low frequency measurements, however, have suggested 
the presence of significant dynamics that appear to be 
strongly correlated with the anomalous structural prop- 
erties^ While it is tempting to attribute these long 
timescale dynamics to the diffuse scattering associated 
with the PNR, several studies using thermal neutron 
spin-echo and cold neutron triple-axis techniques, which 
provide excellent energy resolution, have reported the 
diffuse component to be resolution-limited (and hence 
static) at all temperatures ! 11 ' 28 The supposed static na- 
ture of the diffuse scattering cross section has led to some 
models that rely solely on strain or defects to describe the 
cross section and reconcile it with the dynamics observed 
in infra-red and dielectric measurements. On the other 
hand, other neutron studies have found some evidence 
of low-energy quasielastic scattering that has been sug- 
gested to be correlated with the dielectric properties^r— 
However, because the same low-energy inelastic cross sec- 
tion also contains significant acoustic phonon scattering, 
the basic nature and existence of any quasielastic compo- 
nent remains unclear^ For this reason it is highly desir- 
able to study the diffuse scattering in lead-based relaxors 
with an energy resolution that is comparable to that of 
frequency-dependent dielectric measurements so that a 
direct comparison between the results from these differ- 
ent techniques can be made. 

We have measured the diffuse scattering using cold 
neutron backscattering and spin-echo techniques in an 
effort to characterize the dynamics of the PNR in the re- 
laxor PMN. The large dynamic range offered by both 
techniques allows a direct comparison with dielectric 
measurements, which probe dynamics on timescales less 
than ~ 1 GHz. While our backscattering data repro- 
duce previously published results concerning the static 
nature of the diffuse cross section, our spin-echo exper- 
iments reveal that the diffuse scattering is in fact de- 
scribed by two components - one static and one dy- 
namic. The static component is onset at 400 K, where 
the ferroelectric-active, soft phonon reaches a minimum 
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FIG. 1: a) The square of the soft, transverse optic mode 
frequency is shown as a function of temperature (data taken 
from Refs. Oil and OS. 6) The frequency dependence of the 
dielectric constant of PMN (data taken from Ref. |26| . c) The 
inverse of the 100-kHz dielectric susceptibility as a function of 
temperature. The dashed line represents a high temperature 
fit and the Curie- Weiss temperature (0) is indicated. The 
data were taken from Ref. |l9|. 

frequency. The temperature dependence of the static 
component matches well with our backscattering results 
and with previous data using coarser energy resolution 
performed on triple- axis spectrometers ^ The dynamic 
component is well described by a single relaxational de- 
cay time (r) on the order of ~ 0.1ns. The temperature 
dependence of the relaxational time is described by a 
simple Arrhenius law with an activation temperature of 
HOOK. We suggest that T f ~ 200 K, which is the so- 
called Vogel-Fulcher temperature, is the temperature at 
which the PNR freeze entirely, and that 400 K is corre- 
lated with the onset of static, short-range ferroelectric 
correlations (PNR). The dynamics also reflect the high- 
temperature scale where previous dielectric and optical 
index of refraction measurements have suggested the on- 
set of fluctuating polar nanoregions, often referred to as 
the Burns temperature. 

II. EXPERIMENTAL DETAILS 

Our neutron scattering experiments were conducted on 
the IN10 backscattering and IN11 spin-echo spectrom- 
eters located at the Institut Laue Langevin in Greno- 
ble, France. IN10 consists of a large, vibrating Si(lll) 
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coding the neutron's speed into the Larmor precession of 
its nuclear magnetic moment in a well controlled, exter- 
nally applied magnetic field. I(Q, t) is the spatial Fourier 
transform of the Van Hove self correlation function G(r, t) 
which, essentially, gives the probability of finding a parti- 
cle after time t at a radius r around its original position 34 . 
Furthermore, I(Q,t) is the frequency Fourier transform 
of the scattering function S(Q,lj), which is what is mea- 
sured with neutron backscattering spectroscopy. We used 
the NSE spectrometer IN11 in its high resolution set-up 
"IN11A". The PMN sample is a 3cc crystal grown us- 
ing the Bridgeman technique described elsewhere^ The 
sample has a room temperature lattice constant of 4.04 A 
and was aligned in the (HKO) scattering plane for all mea- 
surements. 



III. RESULTS AND DISCUSSION 
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FIG. 2: a) Diffuse scattering intensity contours measured 
near Q=(100) on the DCS spectrometer located at NIST; the 
circles illustrate where in reciprocal space the measurements 
on IN10 and IN11 were conducted, b) Backscattering spectra 
measured at several temperatures and summed over the Op- 
positions indicated in red. The error bars are related to the 
square root of the total number of neutron counts. 



monochromator that Doppler shifts incident neutrons 
with energies near E{ — 2.08 meV and directs them on to 
the sample. Neutrons scattered from the sample are then 
energy analyzed by a large bank of Si(lll) crystals, which 
backscatter neutrons having energy Et — 2.08 meV 
through the sample and onto a series of detectors. The 
dynamic range for this spectrometer is ±10 /ieV, and the 
elastic energy resolution is SE = 0.5 /ieV (half-width). 
Neutron spin echo (NSE) spectroscopy differs from other 
neutron spectroscopic methods in that it measures the 
real part of the normalized intermediate scattering func- 
tion 5i[/(Q, t)/I(Q, O)}^ 3 -, where Q is the total momen- 
tum transferred to the sample. This is achieved by en- 



We first discuss our measurements of the static diffuse 
scattering cross section, which is associated with short- 
range polar correlations. Fig. [^fr) shows the geometry 
of the diffuse scattering intensity contours in PMN mea- 
sured near Q = (100) at 200 K. These data were obtained 
using the Disk Chopper Spectrometer (DCS) located at 
NIST and have been published and discussed elsewhere*i2. 
The red and white circles indicate the positions in re- 
ciprocal space that we studied using the backscatter- 
ing (1N10) and spin-echo (IN11) techniques, respectively. 
Panel b) shows various constant-Q scans measured using 
IN10. These data have been corrected for a background 
measured at 550 K, above the onset of any diffuse scat- 
tering component. The data in panel b) are the sum of 
the intensities measured at the two Q-positions indicated 
by the red circles labeled IN10 in panel a). As can be 
seen, both points are located far away from the (100) 
Bragg peak position and therefore are not contaminated 
by changes in the Bragg peak intensity as a function of 
temperature. The lines in panel b) are fits to a Gaus- 
sian with a width fixed to the energy resolution of IN10, 
which was measured using the incoherent scattering cross 
section from a Vanadium standard. The data in panel 
b) illustrate the onset of the static (resolution-limited) 
component of the diffuse scattering cross section. These 
data also show the absence of any measurable dynamics 
out to energy transfers of about 10 /ieV between 10 K 
and 550 K (as indicated by the lack of any intensity or 
temperature dependence). However the dynamic range 
probed by backscattering is very limited, and so it is 
desirable to investigate the inelastic properties of PMN 
with a technique covering a broader dynamic range while 
maintaining an excellent energy resolution. Indeed, it is 
possible that the spectral weight gathering within the 
elastic (E — 0) line at low temperatures may result from 
a rapid evolution of the dynamics. To investigate this 
possibility, we have used spin-echo spectroscopy, which 
covers a large dynamic range in time simultaneously. 

The spectra measured using spin-echo (NSE) are plot- 
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FIG. 3: Time dependence of the normalized intermediate 
scattering function $l[I(Q,t)/I(Q,0)] at various tempera- 
tures. Both dynamic and static components are observed for 
200 K < T < 450 K, but only the static component is present 
at 200 K. The error bars are related to the square root of the 
total number of neutron counts. 



ted in Fig. [3] between 200 K and 450 K and characterize 
the intermediate scattering function di[I(Q, t)/I(Q, 0)] as 
a function of t. At 450 K and above, the data are con- 
sistent with a flat line, thus indicating no time depen- 
dent dynamics are present within the range accessible 
with the NSE technique. The error bars are large for the 
450 K data set because the diffuse scattering intensity is 
extremely weak at this temperature. The relatively fea- 
tureless spectrum at 450 K contrasts with that measured 
at 350 K, which clearly shows a systematic decay with 
time, thereby illustrating the presence of a dynamic com- 
ponent. At lower temperatures the decay time increases 
until at 200 K only a flat, ^independent response is ob- 
served, which indicates that the normalized intermediate 
scattering function is purely static within the limits of 
the instrumental resolution, i. e. $t[I (Q , t) / 1 (Q , 0)] = 1. 
The data in Fig. [3] therefore reflects an interplay between 
static and dynamic components of the diffuse scattering 
as the temperature is varied. Based on the direct connec- 
tion between the diffuse scattering and the short-range 
polar correlations observed in PMN, we attribute this in- 
terplay to the presence of both static and dynamic polar 
nanoregions. 

To quantify the static and dynamic components 
we have fit the spectra to a single relaxational form 
n[I(Q,t)/I(Q,Q)] = a + (1 - a)e-^ T . The constant 
a represents the fraction of intensity that is static on 
the timescale of the measurement, and r represents the 
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FIG. 4: Temperature dependence of the static and dynamic 
components extracted from fits to the data in Fig. [3] a) 
compares the static intensities measured with SPINS (5E = 
100/iteV), IN 10 (SE = 0.5 fieV) and IN11. Panel 6) plots the 
dynamic component of the intensity as a function of temper- 
ature. The value of T ma x from Ref. [2^ measured at 1 GHz 
is represented by a vertical line. We emphasize that all data 
was taken on the same PMN sample. 



characteristic decay time of the dynamics. Based on 
the fits shown in Fig. [3J we can extract the fraction of 
the raw uncorrected intensity associated with dynam- 
ics (1 — a) and statics (a). A distribution of decay 
times can be incorporated into the fit by using the form 
$t[I(Q,t)/I(Q,0)} = a + (1 - a)e-('/ T ) . While there 
is physical justification for a distribution of decay times 
associated with the dynamics of the polar nanoregions, 
it is not clear what the value of j3 should be for this 
case. Our choice (J3 — 1) is justified by the fact that our 
data are well described by a single timescale over a broad 
temperature range. We interpret the single decay time 
extracted from the fits to be the average fluctuation time 
of the polar nanoregions. 

The static and dynamic components extracted from 
Fig. [3] have been corrected for the IN11 instrumental 
resolution (7(Q,0)) and plotted as a function of tem- 
perature in Fig. 3J Panel a) compares the elastic (or 
static) intensity measured using three different instru- 
ments and experimental resolutions, and panel b) displays 
the dynamic component normalized to have a maximum 
value of unity. The data taken on SPINS, a cold neu- 
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tron triple-axis spectrometer located at NIST, were mea- 
sured at Q=(0-025, 0.025, 1.05) with an energy resolution 
of SE — 100 fieV. The backscattering (IN10) data were 
measured at the Q indicated in Fig. [2jz) and represent 
the intensity of a Gaussian function of energy fit to the 
elastic peak. The NSE data represent the static parame- 
ter (a) extracted from fits to the NSE spectra described 
above multiplied by I(Q,0). Hence the data plotted in 
Fig. [3}i) equals I(Q, 0) x a. All of the data have been 
normalized to unity at 200 K. The static component de- 
rived from the NSE fits described above agree well with 
the static component derived from backscattering. The 
data do not agree as well with the intensities derived 
from poorer energy resolution measurements on SPINS, 
which used a fixed final energy Ej — 4.5 meV, and are 
suggestive of a dynamic component over this region that 
has been integrated over (i. e. lumped into the elastic 
channel) by virtue of the fact that the SPINS instrument 
provides a substantially poorer energy resolution than do 
the IN10 backscattering and IN11 NSE spectrometers. 

From Fig. @Ji) the onset of the static portion of the dif- 
fuse scattering appears between 400 K and 450 K. This 
is significant when compared to dielectric and other 
neutron inelastic scattering results^ The Curie- Weiss 
temperature derived from high temperature dielectric 
data (Fig. [T]:)) is 400 K 1 — and coincides exactly with the 
temperature at which the ferroelectric-active, soft mode 
reaches its minimum energy (Fig.rijj))^ The 400 K tem- 
perature scale also matches the temperature at which 
strong deviations in the coefficient of thermal expansion 
in PMN are observed^ Based on this result we interpret 
the onset of the static component of the diffuse scatter- 
ing as the freezing of PNR. The short-range nature of 
these frozen regions is evident from the broad nature of 
the diffuse scattering cross section in momentum space, 
which was illustrated in Fig. [2J 

The dynamic component illustrated in panel b) equals 
I(Q,0) x (1 — a). On cooling the dynamic component 
increases, peaks, then decreases; this is consistent with 
dynamics entering and leaving the time window probed 
NSE. The maximum intensity of the dynamic compo- 
nent occurs near ~ 325 K, which is agrees well with the 
temperature T max at which the peak dielectric response 
occurs when measured at a frequency of 1 GHz (see the 
vertical line in Fig. W?))- We emphasize that T max is 
not a meaningful temperature scale because it depends 
on the energy resolution and technique used to measure 
it. Rather, T max is a temperature that is characteris- 
tic of the dynamics and the frequency (or timescale) at 
which they are probed. We also note that contamina- 
tion of our data from acoustic modes, which have plagued 
other studies, is unlikely as these phonons reside at much 
larger energy transfers (or at shorter times ~ 1 THz) than 
are probed here. Also, the acoustic phonon intensity in- 
creases with temperature as required by the Bose factor; 
this does not reflect the trend observed in our experi- 
ment. 

There is little dynamic spectral weight at 200 K where 
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FIG. 5: a) The decay time r is plotted as a function of 
temperature. The solid line is a fit to an Arrhenius law 
(r = T e u/kbT ) with 17=1100 K. b) The fitting parameter 
a is plotted as a function of temperature. 



the normalized NSE spectra (Fig. nearly reaches unity 
and becomes flat. These data can then be understood in 
terms of two temperature scales. The first is an upper 
transition temperature near 400 K where static, short- 
range, ferroelectric correlations are onset. This tempera- 
ture coincides with the Curie- Weiss temperature derived 
from high-temperature dielectric data as well as the tem- 
perature where the ferroelectric-active, soft transverse 
optic mode reaches its minimum frequency. A second 
important temperature scale exists around 200 K. This 
temperature is defined as that where the short-range po- 
lar correlations are truly static. This temperature is most 
clearly reflected in electric field studies, which show that 
the diffuse scattering can be suppressed by an electric 
field only below ~ 200 K. This is also the transition tem- 
perature Tf where dielectric measurements under a pol- 
ing electric field observe a sharp peak characteristic of 
the presence of long-range, ferroelectric correlations. 

We now discuss the temperature dependence of the 
decay time r and the parameter a. Both r and the 
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static component of the normalized intermediate scat- 
tering function 3i[I(Q, t) /I(Q, 0)] are displayed in Fig. [5] 
Panel a) shows the temperature dependence of r, which 
we have chosen to fit to the simple Arrhenius equation 
t = T e u / kbT with U = 1100K±300K. This value for 
U was also confirmed through a universal fit of all the 
data to a single exponential decay following the Arrhe- 
nius equation. Studies of spin-glasses, which exhibit anal- 
ogous slow dynamics, have often described the observed 
relaxation rates in terms of a power law or an exponen- 
tial, Arrhenius-type equatiom 37 ' 38 Regardless of the ex- 
act form, all of these equations describe the dynamics 
in terms of some form of energy barrier and therefore 
have very similar physical interpretations. Some of these 
descriptions for the temperature dependence of the relax- 
ation rate contain a critical temperature as, for example, 
in fits using a Vogel-Fulcher type analysis. More com- 
plicated fits introducing more fitting parameters beyond 
the simple Arrhenius law are not justified by our data. 
The value of the activation energy U = 1100K±300K 
is in reasonable agreement with that (E a ~ 800 K) de- 
rived from infrared and infra-red techniques^ It is also 
in agreement with the temperature at which the lowest- 
energy, transverse, optic mode is observed to become soft 
and broad in energy^ 

Panel b) of Fig. [5] displays the value of a obtained 
from the fits shown in Fig. [3J The value of a varies 
smoothly from 0.5 to 1 with decreasing temperature. We 
do not attribute any significance to the high tempera- 
ture saturation value of 0.5 for two reasons. First, the 
normalization value I(Q, 0) is based on an energy inte- 
gration over the time range of the spectrometer. The 
value of 0.5 could thus be the result of an incomplete en- 
ergy integration (and hence determination of S(Q, 0)) at 
high temperatures. This is corroborated by cold neutron 
spectroscopy data, which provide evidence of the pres- 
ence of dynamics at high temperatures on the timescale 
beyond the resolution of our spectrometer.— Secondly, 
when the measured intensity is taken into account, as 
done in Fig. HJ the results are consistent with no or lit- 
tle spectral weight residing in the static channel at high 
temperatures. The factor of 0.5 does illustrate the pres- 
ence of very long timescale dynamics at high tempera- 
tures. The dynamics are beyond the time window of the 
current experiment, but we speculate that the dynam- 
ics could result from short-range chemical order found in 
Ref. 123 . Further work involving longer timescales are 
required to confirm the origin of this extra component. 
The temperature dependence of a illustrates that all of 
the intensity coming from the diffuse scattering below ~ 
200 K is static. This temperature scale for freezing of the 
polar nanoregions agrees with other techniques. NMR T2 
measurements, which probes fluctuations on the order of 
MHz, have found evidence of an anomaly near this tem- 
perature and have attributed it to freezing of the polar 
correlations^ 

Our results contradict claims based on thermal spin- 
echo measurements presented in Ref . [ill where the diffuse 



scattering is concluded to be solely static. Those mea- 
surements were conducted on a thermal spin-echo ma- 
chine with considerably coarser energy resolution and a 
smaller dynamic range than was used here. We therefore 
believe that these previous measurements only probed 
the static component found in our study. We note that 
the results presented here are consistent with data ob- 
tained from backscattering, in particular the IN10 data 
discussed earlier. Backscattering probes a very narrow 
dynamic range of ~ 1 ns; by comparison, neutron spin 
echo covers several orders of magnitude in time. It is ap- 
parent from Fig. [3] that a much broader dynamic range 
is required to observe both the dynamic and static re- 
sponse. We conclude then that backscattering is not sen- 
sitive to the dynamic component because the dynamic 
range accessed by that technique is too narrow. 

The presence of a dynamic component of the diffuse 
scattering and a favorable comparison with frequency- 
dependent dielectric data strongly link the diffuse scat- 
tering cross section to the presence of dynamic polar 
nanoregions. Our data also clearly illustrate two key 
temperature scales. A high temperature scale, where the 
static component appears (~ 400 K), and a lower tem- 
perature scale (~ 200 K), where no dynamic component 
is observed and only static correlations exist. We do not 
attribute any physical significance to T max , which we 
believe is associated with the dynamics and not to any 
transition temperature. Our interpretation of the up- 
per temperature of 400 K is also in agreement with the 
ideas presented in Ref. [Tt], which further offers an alter- 
native understanding of the Burns temperature at 620 K 
in terms of dynamics. 

Recent reports of a third temperature scale known as 
T* can also be understood in terms of the dynamics. 
Ref. EH describes a new upper transition based on de- 
viations in the thermal expansion coefficient. Such a 
temperature scale may be related to the dynamics and 
could also be associated with the near-surface effect mea- 
sured in single crystal samples of the lead-based relaxors 
PMN and PZN, 4,42 Indeed, previous studies of the coef- 
ficient of thermal expansion in PMN using neutrons have 
found a strong variation as a function of depth. While 
these speculations are inconclusive, further study is re- 
quired to understand the origin of these upper transition 
scales. Dielectric and infra-red work in Ref. describe 
a third temperature scale T*, which is 400 K, and we 
interpret this temperature as the onset of static polar 
correlations. Ref. |43| also descrobe a high temperature 
scale near 600 K. This could be related to the dynamics 
as evidenced by the large activation energy derived in our 
current experiment. 

The interpretation of the dielectric and neutron data 
in terms of two temperature scales is in agreement with 
random field models previously proposed to explain the 
relaxor ferroelectric phase diagram. References j44j-(46| in- 
terpret relaxors in terms of a 3D Heisenberg model with 
cubic anisotropy in a random field imposed by the dis- 
order on the perovskite B site. In this picture, the true 
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transition temperature would correspond to 400 K where 
the soft, transverse optic phonon mode energy reaches 
a minimum. However, this model also predicts a sec- 
ond (lower) temperature scale when the energy scale im- 
posed by the cubic anisotropy becomes important. It was 
suggested that only below this temperature could long- 
range, ferroelectric correlations develop in the presence of 
a random field. This model is in broad agreement with 
the data presented here, which can be described using 
just two temperature scales. 

We have presented a high resolution investigation of 
the diffuse scattering in PMN as a function of temper- 
ature that sheds new light on the underlying behavior 
of the polar nanoregions. We have proven the existence 



of a dynamic component to the diffuse scattering and 
that the temperature dependence is described by an in- 
terplay between dynamic and static components. We find 
that the static component appears between 400 K and 
450 K, which is in excellent agreement with the Curie- 
Weiss temperature derived from dielectric data and also 
the minimum in the soft mode energy. From our spin- 
echo spectra, we have also derived a dynamic fraction and 
shown that it peaks near the temperature T max found in 
frequency-dependent dielectric constant measurements. 
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